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Abstract 

In three-dimensional magnetic configurations for a plasma in which no closed field line or magnetic null exists, no 
^ magnetic reconnection can occur, by the strictest definition of reconnection. A finitely long pinch with line-tied 

boundary conditions, in which all the magnetic field lines start at one end of the system and proceed to the opposite 
end, is an example of such a system. Nevertheless, for a long system of this type, the physical behavior in resistive 
QQ magnetohydrodynamics (MHD) essentially involves reconnection. This has been explained in terms comparing the 

i-H geometric and tearing widths yji2j|. The concept of a quasi- separatrix layerfT,*?! was developed for such systems. 

In this paper we study a model for a line-tied system in which the corresponding periodic system has an unstable 

^ f-H tearing mode. We analyze this system in terms of two magnetic field line diagnostics, the squashing factor^\^ and the 

electrostatic potential difference used in kinematic reconnection studies ['S', 7 |. We discuss the physical and geometric 
significance of these two diagnostics and compare them in the context of discerning tearing-like behavior in line-tied 
modes. 
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'c/3 1. Introduction 

>> 

4^ It has been argued that in three dimensions hyperbolic closed field lines (X-lines) f6l or magnetic nulls ['6^, "81 are 

^ Oh required for magnetic reconnection. In two dimensions, reconnection has also typically been discussed in terms 
of nulls (zero guide field) or X-lines (nonzero guide field). See Ref. [9i| and references therein. In this strictest 
i-H definition of reconnection, the argument is this: in ideal magnetohydrodynamics (MHD) these structures can lead to 

^ singularities which are then resolved by non-ideal effects, and in the absence of such structures, no singularities occur. 

In the presence of small electrical resistivity r] , these singularities become current sheets of finite but small thickness 
dependent on r] . For our purposes here, we will describe the non-ideal effects as arising from finite electrical resistivity 
r], i.e. we consider reconnection in resistive MHD[^In this strict definition, the question of whether reconnection occurs 
is in the context of the limit r] ^ 0, or the Lundquist number S ^ oo. 

In a line-tied system of finite length — L <z<L with fields having > 0, all field lines enter at z = —L and exit 
^-H at z = L. Therefore, closed field lines cannot form. (We will speak only in terms of closed field lines; in order for 

i-H nulls to form in this geometry, all components of B would have to change sign in the region.) So in the strict sense 

^ described above, no magnetic reconnection can occur. However, in Refs. [1 , 2J line-tied modes which are tearing modes 

• in cylindrical geometry [for equilibria with BQ(){r), Bzo{r)] with L = oo were studied. In these L = oo modes, behaving 

as e^^^^^^^, there is a mode rational surface consisting of closed field lines in the equilibrium, and the perturbed fields 
^ have hyperbolic and elliptic closed field lines (islands). The conclusion of these papers was that the modes still behave 

as tearing modes (with growth rates 7 proportional to a fractional power of T]) in a sufficiently long line-tied system 
L > Lcrit- The condition L > Lent is defined as the range over which the tearing width Wt of the mode is greater than 
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^However, there is no fundamental reason why these concepts cannot apply to reconnection with other non-ideal effects, such as electron inertia 
and hyper-resistivity. 
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the geometric width w^, the width of the mode associated with the hne-tying boundary condition [Wg oc 1/L]. For 
L < Lcrit, when Wt <Wg, the modes no longer resemble tearing modes; they have no recognizable tearing layer and 
ycx: 7], i.e. they involve global resistive diffusion rather than reconnection. 

These results suggest that in many cases it may make more sense to define reconnection for fixed but large S, rather 
than in the limit S ^ oo. This is particularly true for numerical simulations with limited values of the Lundquist number 
S and for laboratory experiments with relatively low temperatures. Three such experiments suited for reconnection 
studies with line-tying are RWM lfTOl , RSX lfTTIl . and LAPD |[T2l . This point of view is consistent with that expressed in 
Ref. 13] Hi, in which the concept of a quasi- separatrix layer (QSL) was introduced for fields without hyperbolic closed 
field lines (or nulls) as a resolution of controversies associated with generalized magnetic reconnectiorj^ 13^ J4J- In 
fact, the concept of a QSL was used to analyze LAPD in ifTSl . For such magnetic fields (with large but fixed S), the 
geometric aspects of the field lines can act to separate the field lines in a manner qualitatively similar to fields with 
hyperbolic lines. If this separation occurs in a thin enough region, the physics is basically identical to the behavior in 
the presence of a hyperbolic line. However, note that while the condition Wg < Wt may be satisfied for such laboratory 
systems and for simulations with relatively small S, for solar coronal or astrophysical applications this requirement 
may hold only for unrealistically long plasmas. 

In Refs. |[3]-[5l [T6l - [T9i , it was suggested that the most effective measure of a QSL is the squashing factor, which 
measures the squeezing and stretching of field lines. The method of slip- squashing factors 1 19] also takes resistivity (or 
other non-ideal effects) into account. It involves computing the squashing factor plus taking into account the field line 
slippage due to non-ideal effects. For fields which are determined by boundary motions, this slippage is measured by 
comparing the initial field line mapping with the field line mapping at a later time. 

In Sec.|2]of this paper we discuss the squashing factor Q and a second diagnostic, related to the potential 0^ used in 
Refs. in [71. The latter function is the scalar potential required to give E • B = 0, a consequence of Ohm's law in ideal 
MHD. We first compute Q and 0; for two examples of magnetic fields, namely a doublet-like field and a field given by 
an equilibrium BQo{r), Bzo{r) plus a single tearing mode, both in a finite region — L < z < L. We also compute 0^, the 
scalar potential required to satisfy E • B = r]j • B in resistive MHD. In Sec. 3 we review the treatment of line tied modes 
in the two-mode approximation\20\, which applies for long plasmas. In Sec. 4 we study QSLs in a model representing 
growing tearing modes with line-tied boundary conditions in the two-mode approximation. We show the squashing 
factor Q as well as the potentials 0^ and 0^ for modes of various amplitudes, in the two cases Wt > Wg and Wt <Wg. 
In Sec. 5 we summarize and discuss our results. In the Appendix we discuss other mathematical issues related to the 
squashing factor. 

2. Reconnection diagnostics 

In this section we describe two diagnostics that can be used to identify reconnection and QSLs in systems with 
finite resistivity. The first diagnostic is the squashing factor, defined in Refs. HSlfTBlTlS L and the second is the potential 
difference which is required to satisfy either the ideal MHD relation E • B = (see Ref. |6|) or the resistive relation 
E • B = r7j • B. Each of these diagnostics require computing the magnetic fields lines which run the length of the system, 
from z = —L to z = L. 

2.1. Squashing Factor 

The first diagnostic we consider is the squashing factor. Integrating field lines from one end of the system to the 
other defines a coordinate mapping M : x ^ X, where x is the starting point of the line in the z = —L plane, and X 
is the ending point in the z^L plane. (In the more general case, we map from the set with normal field 5„ < to 
the set with 5^ > 0.) Geometrically, stretching and squashing of flux tubes indicates a potential for reconnection, 
and the Jacobian / of the mapping M gives us information about these processes. The overall expansion of flux 
tubes is not associated with reconnection, but rather is related to variation of the guide field strength B^. Specifically, 
flux conservation implies Bz{x\.,X2-, —L)dx\dx2 = Bz{X\.,X2^L)dX\dX2, so that the Jacobian determinant D = det(/) 
satisfies D = 5^(xi ,X2, —L) /5^(Xi ,X2,L). The mapping M takes a flux tube with a circular cross-section and "squashes" 



^The controversies involving GMR related to whether localized current had to be from reconnection or could be due to double layers, Pfirsch- 
Schluter currents, or other sources. 
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it so that its cross-section is elliptical. In order to measure the stretching and squashing while compensating for the 
expansion of flux tubes if const., a natural quantity to consider |[3]-[5l [T6HT9ll is the aspect ratio of the elliptical 
cross section of the flux tube at z = ^0 




Figure 1: Squashing of flux tube cross-section by the mapping M. The major and minor axis of the eUipse have lengths pmax = max(|5X|) and 
Pmm = min(|5X|) (where p/ are the singular values of the matrix /), and so the aspect ratio of the ellipse is = pmax/pmin' This ratio removes the 
effect of variations in B^, which would only give an overall expansion or contraction of the flux tube. 

Consider two field lines whose initial points are separated by the tangent vector 5x. Their endpoints will be 
separated by 5X, and the length of 5X compared to the length of 5x will tell us about the stretching of a flux tube 
containing the two field lines. The lengths of these two vectors can be directly compared since the Jacobian matrix / 
(Jij = dXi/ dxj) maps tangent vectors, i.e. 5X = /5x. Thus, their lengths are related by 

(5X, 5X) = (/5x,/5x) = (5x,/^/5x). (1) 

The symmetric positive definite matrix J is the covariant metric tensor on the surface at z = —L derived from 
a Euclidean metric tensor on the surface dX z = and its eigenvalues A±, obtained from the Rayleigh quotient 
(5x,/^/5x)/(5x, 5x), tell us how the shape of a flux tube changes]^ In fact, the ratio of the singular values L21J 
P± = \/^ equals the aspect ratio of the ellipse. See Fig. Ill This ratio is easily shown to obey 



R=^ = W+ \ (?) -1> where Q = T/D = Tr{J^J)/\dct{J)\. (2) 

Pmin ^ V V / 

This quantity 2, which also equals 7? + 1/7?, is called the asymptotic squashing factor i 1 8 J . or just squashing factor. In 

regions where Q is large, R = Q—\/Q^ ^ Q- It was argued in Refs. lISHSl [T6HT9ll that these regions, where the 

stretching and squashing are large, are candidates for the occurrence of reconnection. Such regions of large Q are 
called quasi- separatrix layers. We are interested in only these regions and, following Refs. |[3ti5l [T6HT9 K we will thus 
focus on Q rather than R. In the next section, we give an example which illustrates how Q becomes larger and more 
concentrated for larger L. 

2.2. Squashing Factor Example: Doublet Field 

2 2 4 

Consider the vector potential with A^^ = —BqX2, A^^ = 0, and = ^ — ^ + ^ . This defines the magnetic field 

Bjc^ X2, = ^1 - , B,=Bo. (3) 

The equations for field lines are then 

dxi/dz = X2/Bo, dx2/dz = {xi -x\)/Bq. (4) 

These equations are in canonical form, with Hamiltonian H = A^(xi ,^2) /Bq and D = I (since Bz= Bo = const.). The 
contours of A^ for this field are shown in Fig.|2^. Integrating the field lines from z = —L to z = L with a symplectic 
integrator determines an exactly area preserving mapping (xi,X2,— L) (Xi,X2,L). The Jacobian matrix of this 



^In the Appendix, we discuss the possibility of transforming (xi ,X2) and (Xi ,^2) to canonical variables, for which D = 1. 
^The quantity Tr(/^/) was described as the norm of the displacement gradient tensor in | 3 |. It can also be written as \VXi p + IVX2P and is 
actually the square of the Frobenius norm of /. 
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mapping is estimated numerically by tracing field lines for a grid of initial points x, and taking the difference of the 
ending points X. Since this model is two dimensional, we can change the length L without modifying the model, as we 
do in the 3D models in Sec.js] The squashing factor 2 as a function of (xi ,^2) is shown in Fig.[2]for lengths L = 1, 2, 
and 2.5. Note that for L = 2, 2 is peaked in an elliptically shaped region near the part of the stable manifold near the 
X-line. As L increases, Q becomes much more peaked and concentrated near the stable manifold and further along 
the stable manifold. The squashing factor Q is shown in Fig. [2]: for L = 2.5 as a function of the variables at z = L, 
namely (Xi,X2). The large values of Q in this figure are concentrated on the unstable manifold because this same 
quantity Q{X\^X2) could also have been computed by initializing at (Xi ,X2,L) diiz = L and integrating backwards to 
(xi,X2,-L). 

In Fig.|3]we show a three dimensional plot of the contour where Q = 1000 for L = 2. The concentration along 
the stable manifold for z = —L and along the unstable manifold for z = L is evident. This surface has the form of a 
hyperbolic flux tube, which arises in models of reconnection in solar physics, and the study of QSLs. |[T7l[22ll23l 




Figure 2: Squashing factor ,^2) for the model of Eq. |3j, with 5o = 1 and (a) L = 1, (b) L = 2. In (c), the squashing factor for L = 2.5 is shown 
plotted as a function of ending point X. In (a), contours of are shown, with arrows indicating the stable and unstable manifolds of the X-line at the 
origin. In this 2D case, the stable and unstable manifolds coincide to form a separatrix. The quasi-separatrix layer traces out the stable manifold 
when plotted as a function of initial point x [(a) and (b)] and the unstable manifold when plotted as a function of the ending point X in (c). 




Figure 3: 3-dimensional visualization of the quasi-separatrix layer for the model of Eq. ([Sj. The surface is formed by the field lines which have 
e= 1000forL = 2. 



2.3. Scalar Potential Difference A(j) 

The second diagnostic we consider in this paper is the potential difference between the two ends of the field lines. 
We calculate A(j) = A(j)i from the ideal MHD Ohm's law E + vxB = Oto give zero parallel electric field £"11 oc E • B, 
and calculate A(j)r similarly from the resistive MHD Ohm's law E + v x B = r]j. (Again, we consider resistive MHD 
just for concreteness.) The first approach, for ideal MHD, was used to analyze three dimensional reconnection in the 
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presence of hyperbolic closed field lines or nulls in Ref . (S] . In this paper, we compute both A^^ and A(j)r by taking the 
magnetic field B(x,^) and the inductive electric field —dA/dt computed from resistive MHD. We solve E • B = and 
E • B = r7j • B along field lines, namely 



B-V0, = -B— ; B-V0, = -B — -r7j-B, (5) 

for the scalar potentials 0/, 0^ and, for both, evaluate the difference between the potential at the two ends of the field 
line at z = —L and z = L. The potential 0^ is exactly the single- valued scalar potential related to the fields obtained 
in resistive MHD (in a specific gauge) and represents the electrostatic potential due to small (quasineutral) charge 
variations required to balance the inductive electric fields of the mode locally along the field lines, in the presence of 
the resistive term T]] • B. 

In computing 0^, we analyze the fields generated by resistive MHD in the context of ideal MHD, in order to 
determine when ^ In most regions of the plasma, the potential 0^ is also a single- valued function, again due to 
movement of electrons along B to cancel the inductive electric field. In regions where rfyBis small, A(j)r ~ A0/. On 
the other hand, if there is a closed field line in the system and if dA/dt integrates to a nonzero value along this line, 
then mulit-valuedness of 0^ will occur at points on the line and singular behavior will occur on points that go to the line 
(the stable manifold of the hyperbolic X-line), because the associated flux change through any surface bounded by the 
line cannot be cancelled by a single- valued electrostatic potential]^ In the line-tying geometry of this paper, we do not 
expect to encounter true singularities in 0^. However, we do find that 0;, like the squashing factor Q, can be large and 
strongly localized, which is indicative of reconnection-like behavior, e.g. current sheets. This kind of behavior, with 
peaked but nonsingular values of 0/, was seen and discussed in |7|. 

It is useful to compute both (j)i and (^r, since the difference A^^ — A(j)r is simply the integral of r]] • B oc r] j|| = £"11 
along the field line. Parallel currents are associated with reconnection effects, and this method allows us to not only 
compute the parallel currents (through A^/ — A0^), but it also gives us quantities (A^^ and A(j)r) with which to compare 
the parallel currents. If A(j)r ~ A(j)i, then the effect of the parallel currents are small compared to the inductive electric 
fields. In regions where A^^ — A(j)r is large compared to A0/, however, the currents - and thus reconnection effects - are 
important. The difference A(j)i—A(j)r is equal to the quasi potential S of Ref. |24|. However, we find that comparing the 
two separate potentials A^^ and A0^ is more informative than considering only S = A^^ — A(j)r. 

As diagnostics for reconnection, there is one fundamental difference between the squashing factor Q and the 
potentials 0/, (j)/. The former involves only the structure of the magnetic field at a specific time, whereas the potentials 
(^/, (l)r obtained from Eq. ^ involve the structure of the magnetic field and the inductive electric field —dA/dt. That is, 
an area of large and localized Q only indicates the possibility of having reconnection, whereas a large and localized 
value of A(j)i shows flux changes and field line topology indicative of reconnection. 

2.4. Ideal MHD Scalar Potential Example: Single tearing Mode 
We use the compressible zero-j3 visco-resistive MHD equations: 



dtP = -V-(pv) (6) 

p(9^v + pv-Vv = jxB + vpV^v (7) 

E + vxB = r]j (8) 

VxE = -dtB (9) 

VxB = j (10) 

V-B = 0, (11) 



in a cylinder with — L <z<L and < r < r^^, with r^^ = 2. 

The equilibrium used was the equilibrium of 1 1 1, for which the fastest growing mode in an infinite cylinder is a 
tearing mode rather than an ideal MHD mode. This force-free equilibrium is specified by its axial current density 



^Note that a gauge change A ^ A + V;^, ^ ^ — dx/dt with x continuous does not affect the singular nature of as found from Eq. js}. 
Similarly, in integrating Eq. |5| from z = —L to z = L, a smooth initial value (p{r,Q^z^ —L) is irrelevant. 



Ao(^) = 2/(1 + r'^) with K* = 6, and 5^o(^) determined by force balance with an integration constant 5^o(^ = 0) = 5, 
with Po = 1 and vq = 0. This equihbrium is "tokamak-hke" in the sense of having an increasing profile of the field line 
pitch ii{r) = rBzo{r)/Beo{r) and^^o > Bqo- 
The linearized MHD equations are 



dtY = [V X B - A (r)B] X Bo + v V^v 
dfB = Vx (vxBo-r7VxB) , 



(12) 
(13) 



with B = B(r)^^^^+^^^ and v = v(r)^^^^+^^^. These equations are solvecQas an eigenvalue problem for the growth rate 
7, with m = 1, and k = —0.13507, giving a mode rational surface (where k • Bq = 0) at = 1.2576. For computing 
the field line diagnostics, we have superimposed the equilibrium fields with the perturbed fields multiplied by a small 
amplitude The system length L was taken to be L = 2071 /k. This length and the value of k were chosen for 
comparison with the results from the two-mode approximation in the next section. 

We have computed A^^ by integrating Eq. ^ along field lines of the field Bq + aiB from z = —L to z = L, using 
a=lx 10~^ (see Fig.|4]). This quantity is concentrated around the mode rational surface at = 1.2576, near both the 
X-line (6 ^ 7r/2), and the elliptic (O) line (0 « 371 /2). On both the hyperbolic and the elliptic closed field lines, \A^i\ 
is large because the perturbation is constant along these lines. For larger aniplitude a = 0.01, the island opens up, and 
\A^i\ becomes localized around the stable manifold of the X-line (see Figs.|5p,c). We have also shown the surfaces of 
the helical flux ;^ (r, + fe) = Xo{r) X{^^ ^ kz) = mA^ — krAg at z = —L. Surface of section points for periodicity 
Lp = 27i/k lie on surfaces of constant X- 
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Figure 4: (a) A(j)i{r, 6), and (b) A0^(r, 6), computed by integrating Eq. js} from z = —L to z = L. Mode amplitude a = 1 x 10~^. (c) Difference 
between A^, and A^^, which is the cumulative parallel current integrated along the field lines. Contours are of the helical flux, which for one mode is 
identical to the surface of section. The white arrows indicate the direction of the stable manifold of the X-line near (r, 0) = (1.25, 7r/2). 



The potential A^^ can also be computed analytically. For infinitesimal 0^ and A we use dO/dz= 1/m(^) the 
equilibrium fields to solve Bq • V0; = — 7B0 • A. With the normalization chosen so that A0{r) and A^(r) are pure 
imaginary, the potential difference iF] 



= 2L 



yBp-A 

yBp-A 

^zO 

yBp-A 

BzO 



dz sin{m6{z) ^kz) 
cos(m0p + m{z ^L)/ii^kz) 



k^m/jl 



sin(m(0p+L/;U)) 



sin((^ + m//i)L) 
{k^m/il)L 



(14) 
(15) 
(16) 



^As in Ref. (25), a small amount of divergence cleaning is used. 

^For even relatively small values of the amplitude a, nonlinear effects could be important, invalidating the assumption that we can superimpose 
the equilibrium and perturbed fields obtained by linear theory. For our purposes in this paper, we will keep a small so that such errors are small. 
^The r dependence of A, Bq, and /i is suppressed here for clarity. 
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where 60 = 3iz = —L. The quantity (m//i(r) = k • Bq/^^o. so A^^ is proportional to a sine function with a peak 
at the mode rational surface r = r^. The width of in r is 



mLii\rs) 



(17) 



Note that the cosine factor in A^^ depends on both and r [through /i(r)], so its zero contours will be tilted in the 
(r, Qq) plane. For this calculation with infinitesimal mode amplitude, there is no information about the hyperbolic and 
elliptic lines, and in fact the island is assumed to have zero width. The width Wm represents the decorrelation due to the 
shear in the magnetic field at the mode rational surface, which is proportional to /i'(r^). This width is consistent with 
the observed radial width of A(^/ in Fig.|4^. Note that in Fig.|4^, the similar prefactor (B • A)/5^ goes rapidly to zero 
outside the mode rational surface, so the amplitude of A0/ falls off more rapidly than might be expected from the sine 
function alone. A plot of A0/ as given by Eq. ( 16) [including theprefactor (Bq • A) /^^ol is nearly indistinguishable 
from the plot in Fig.|4^, which was computed by integrating Eq. ([s^) numerically along field lines of Bq + aB. 
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Figure 5: (a) Ideal and resistive forms of A(|) as in F ig.[4^ ,b, but along the line = 7r/2. Note the large peak in A(|)/ near the mode rational surface, 
(b) A(|)/(r, Q) and (c) A0/(r, Q) — A(j)r{r, 0), as in Figs.|4^ and|4]^, but with mode amplitude a = 0.01. Notice how the structure becomes concentrated 
near the stable manifold of the X-line. 



Fig.|4j5 shows A(j)r, computed from Eq. The oscillations below r =1.2 are the same as in Fig.|4ji, but the 
larger peaks near r = have disappeared, because — 7B • A balances r]} ■ B near r = rs. The difference A(j)i — A(j)r is 
plotted in Fig. [4]:, showing only the voltage drop due to the resistive term near r = rs. In Fig. [5^ we show A^^ and A0^ 
as a function of r for = 7r/2. The large peak in near r = is evident in A^^ but not in A0^. We show A^^ in Fig.[5|), 
as in Fig. [4^ but with mode amplitude a = 0.01. In Fig. [5]: we show the difference A^^ — A(j)r as in Fig. [4]:. We see that 
the minimum and maximum values of A(j)i, at the X-line and the O-line, respectively, are equal and opposite, but the 
values near the X-line are more concentrated along the stable manifold of the X-line. In Fig.[6j we show a plot as in 
Fig. [5^ for L = 252. There is very little difference between A0; and A(j)r for this case. As we will discuss in the next 
section, this is consistent with Wm being larger than the tearing mode width Wt. 



3. Line-tied tearing modes 

For analyzing modes with line-tying, we use the method of Ref. We expand the perturbed velocity and 
magnetic field components ^ as sums over the spectrum: 

ar,e,z,t) = e'"''+yl^a,^,{r)e''\ (18) 

k 



These radial eigenf unctions satisfy the linearized equations (12) and (13 ) and the wavenumbers k and the expansion 
coefficients a^ are chosen so that the solution satisfies the line- tied boundary conditions at z = ±L. If L 1, then 
keeping only two terms k\ , ^2 in the sum gives a good approximation to the line-tied solution. This is the two mode 
approximation. It is valid for large L because — k\ is small for large L and therefore (^) — (^)- See Ref. L20il . 
The calculation proceeds as follows. 
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Figure 6: Ideal and resistive forms of as in Fig.[5j, but with parameters v, T\,ki, and L as in Case II (see Table^. While the difference between 
and A(|)^ is largest near the mode rational surface, it is still small compared to A0/. This is an indication that mode width Wm (Eg. [17) is larger 
than the tearing width Wt for this mode, and thus we expect that the effects of the Hne tying will be important. 



The line-tied boundary conditions at z = ±L correspond to having the tangential components of E and the normal 
component of B be zero at z = -^L. (We describe these modes as line-tied, in spite of the fact that the tangential velocity 
at z = -^L is nonzero due to finite resistivity.) We proceed by treating y as fixed for now, and choose k\ and such 
that in the infinite cylinder results 7(^1) = 7(^2)- Then the line-tying condition leads to a\e^^^^ ^ a2e^^^^ = at z = 
This implies that k2{y) — k\ (7) =l7t/L, giving a dispersion relation for 7. Here / is an integer, but we can take / = 1 
without loss of generality, and for / = 1 we have a\= a2 = a. The length of our line-tied system is then 2L, where 
L = 7t/{k2-ki). 

For this line-tied mode satisfying the two-mode approximation, we calculate the geometric width Wg in the following 
manner|l |. We find the two mode rational surfaces (k • Bq = 0) corresponding to k\ and ^2 by m/ il{ri) -\- ki = 
m//i(r2) +^2 = 0- The geometric width is the distance between these mode rational surfaces: 



: r2 - ri 



{k2-k,)^{n)^ 



mil'{ri) 



(19) 



In Ref. |T| it was shown that this definition is still qualitatively correct even for smaller lengths L, for which the 
two-mode approximation is not valid. It is interesting to note that the single mode width Wm in Eq. p7| equals Wg if the 
line-tying relation k2 — k\ = tt/L is substituted. We shall return to this point in the next section. 

For a large enough, the islands corresponding to the two modes will overlap and magnetic field line chaos will 
occur|^ In order to be able to construct a surface of section plot with fields having a common period in z, we do the 
following: We first assume that ^2/^1 = ^2/^1 is a rational number. In this case the two modes then have a common 
length Lp over which they are periodic, Lp = n\L\ = ^2^2, where Li = 27t/\ki\. Also, it can be seen from the line-tying 
relation above that L = ^Lp/\n2 —ni\. So, if we take |^2 — ^1 1 = 1, then Lp = 2L. From this we conclude that the fields 
are actually periodic over the length of the line-tied system z = —L L. For this special case, the map from z = —L to 
z = L is a surface of section, and we can plot its return map to assess chaos in the periodic system. 

Given the perturbed fields B from this two-mode calculation, we can find A in the gauge where = 0. For the ki 
component, we obtain 



A,i(r) = - rBei{r')dr', 
Jo 



(20) 



~ jfi ~ I ~ 

Aei(r) = — A,i(r) + -B,iW, (21) 



and similarly for the ^2 component. 



^As in the last section, we will allow superposition of states with relatively large a even when nonlinear effects are important. For sure, having 
field line chaos is such a nonlinear effect. 
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The perturbed vector potential for the two mode approximation is thus 

A(r,e,z,0 =^^^^'"^ (ki{r)e'^^'^k2{r)e'^^') . (22) 

The inductive electric field E/ = — is then: 

E/ = -yae^'^'"^^ (^ki{r)e'^'' ^ k2{r)e'^^') . (23) 

This is the inductive field used for computing the scalar potentials A^^ and A0^ [c.f. Eq. ([5])]. This mode behaves in z as 
^i^o^cos(5/:z), where is (^i +^2)/2 and 5^ is k2 — ko. Therefore at r = ro ^ (ri + r2)/2 the mode is nearly constant 
along the field lines for most of the length of the plasma. 

4. Reconnection diagnostics for line tied modes 

Given the two-mode approximation for line-tied modes, we use the perturbed fields B to compute the squashing 
factor Q and the scalar potentials 0/ and 0^ of Eq. ([5]). Two cases are compared: first, the plasma parameters are 
chosen so that the two-mode approximation is fairly accurate, and the tearing width Wt (estimated from the width of the 
perturbed current near the mode rational surface) is larger than the geometric width Wg. The second case has these 
limits reversed, with Wt < w^; in this case, the two-mode approximation is not as accurate, but can be used as a simple 
model for understanding qualitative aspects of the line- tied mode. 



Table 1 : Parameters for the two line-tied cases 

77 V ^2/^1 h k2 7 Wg Wt Li L2 L 

Case I: 19/20 -0.135 -0.128 0.0233 0.047 02 46^5 49^0 465 

Case II: 10"^ 10"^ 11/12 -0.149 -0.137 0.0110 0.084 0.05 42.1 45.9 252 



4.1. Case I: Wt > Wg 

The plasma parameters for this case are t] = 1 x 10~^ and v = 1 x 10~^. The two modes have ^2/^1 = 19/20, 
L = 465, and the values of ki, 7, and Li are given in Table [T] Two different mode amplitudes are shown for this case. In 
Fig.|7j the mode amplitude is a = 1 x 10~^. The surface of section shows two sets of islands with very thin secondary 
resonances between them. The A0; calculation (Fig.|7^) shows structure near the islands for the two modes. The fact 
noted in the last section, just after Eq. ( 19 ), namely that Wm = Wg, implies that the regions of large |A0; | from the two 



island chains overlap. This is true because Wg, as defined in Eq. ( [T9| ), is the radial displacement between the island 
chains. The values of \A(j)i\ are peaked near the secondary resonances with r = ro, between r = ri and r = r2, and 
aligned with the stable manifolds of the X-lines of the two island chains. On the other hand, A0^ (Fig.[7]3) is much 
smaller in magnitude and not localized. The squashing factor is also localized near the islands (Fig.|7|:). Figure [8] shows 
results for the same two modes, but with amplitude a = 5 x 10~^. Note that the surface of section now shows a chaotic 
region linking the islands, and A^^ (Fig.[8|i) and squashing factor Q (Fig. [8]:) are both now more localized in the chaotic 
region and aligned with the stable manifolds of the X-lines. The resistive potential A0^ is smaller in magnitude and 
shows little variation in the island region (Fig. [8})). 

4.2. Case II: Wt < Wg 

The plasma parameters for this case are t] = 1 x 10~^ and v = 1 x 10~^. The two modes have ^2/^1 = 11/12, 
L = 252, and the values of ki, 7, and Li are given in Table [T] For this shorter system, A(j)i and A(j)r are nearly equal 
throughout the plasma (Fig. [9^), meaning that the resistive voltage drop, proportional to A^^ — A0^, is very small. This 
is an indication that the mode has lost its tearing character because of the shorter length of the plasma. This effect 
is quantified in Ref. 1 1 1, where it was shown that the scaling of the growth rate is tearing-like (with 70c a fractional 
power of T]) when Wt > Wg as in the previous case, but the mode exhibits resistive diffusion (with rj) whenever 
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n/2 71 37i/2 In ' 7i/2 71 37r/2 Iti ' 7i/2 71 37r/2 27i 

e e 



Figure 7: Reconnection diagnostics for the two mode approximation with mode amplitude « = 1 x 10~^. All quantities computed with field lines 
mapped from z = —L L. (a) A(j)i(r, 6),(b) A0^(r, 6), and (c) squashing factor Q(r,6). Black dots show the surface of section map, which shows 
two sets of islands, and a secondary resonance. Stable manifolds of the X-lines are indicated by the white arrows. 




Figure 8: Reconnection diagnostics for the two mode approximation with mode amplitude a — 5 x 10~^. All quantities computed with field lines 
mapped from z = —L L. (a) A(j)i{r, 6), (b) A(j)r(r, 6), and (c) squashing factor Q{r, 6). Black dots show the surface of section map, showing field 
line chaos due to the overlap of the two sets of islands. The two X-lines are still near = 7r/2, r = ri and 6 = 37r/2, r = r2, respectively. 
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Wt < Wg. Even in the region where the current is strong (the tearing layer) the finite length is a more important factor 
than the influence of the current. On the other hand, the squashing factor (Fig.[9]3) is quite localized in the region of the 
two island chains, indicating the stretching character of the X-lines of the islands. In this case, the variation of Q in 
this region suggests the potential for reconnection, whereas the potentials A^^ and show that for this case there is 
essentially no reconnection going on. In contrast with 2, the slip-squashing factor of Ref. |[T9l measures the actual 
field line slippage due to non-ideal effects as well as stretching. 




0.5 1 1.5 2 1 1.2 1.4 



r r 

Figure 9: Reconnection diagnostics for the case Wt <Wg. (a) A(|)/, A0^, and a measure of the parallel current integrated along field lines, given by 
— A(j)r (a = 10~^, = 7r/2). In this case, the finite length of the system is more important than the parallel current, as can be seen in (a), since 
A(j)i and A(j)r are nearly the same. In (b) is shown log^olG) for various mode amplitudes a. The peaks in Q are quite evident for large a but not 
noticeable for small a. The importance of the finite length (wt < w^) is not obvious from Q. 



5. Summary and conclusions 

We have constructed a cylindrical model for fields with a growing mode in the presence of line-tying at z = ±L. 
The fields in this model consist of an equilibrium plus perturbed fields E, B of an arbitrary but small amplitude. The 
perturbed fields are computed using the two-mode approximation. The two infinite cylinder modes, with distinct 
axial mode numbers ki, fe, are both tearing modes, so the line-tied mode formed from them should have tearing 
character HJ 121 (i.e. effectively involve reconnection) for large L. Also, it is in this limit that the two-mode approximation 
is valid. Refs. |[l]|2l showed that for shorter plasmas, the modes no longer have tearing character. We have traced field 
lines and computed the squashing factor Q and the electrostatic potentials A(j)i, A(j)r obtained from Eq. ([5]), where the 
inductive electric field is given in terms of the linear modes by E/ = — 7A and yis the linear growth rate. This model, 
with an equilibrium plus linear modes of arbitrary amplitude, is accurate for very small mode amplitudes but clearly 
inaccurate for large amplitudes, for which the nonlinear terms are important. Quantitative results for large amplitude 
nonlinear states require nonlinear resistive MHD simulations, which are outside the scope of the present investigation. 

We discussed the point that the squashing factor or squashing ratio depends on only the geometry of the magnetic 
fields, whereas the electrostatic potentials A(j)i and A(j)r involve the magnetic field geometry plus the inductive electric 
fields. So, large and localized values of Q indicate a magnetic structure with the potential for reconnection; large and 
localized values of A0/ relative to A(j)r indicate electric and magnetic fields for which reconnection should be occurring. 
A comparison between our methods using A^^ and A(j)r and the use of the slip- squashing factor of Ref. LI 91 should be 
interesting, but is also outside the scope of the present study. 

We found very different behavior for A^^-, A0^ in Case I, the long plasma case (Wg < Wt) and Case II, the short 
plasma case (Wg > Wt). In Case I, A^^ is much more peaked in the tearing layer region than A0^, indicative that 
reconnection plays a large role in the line-tied mode in this range. In Case II, A^^ and A(j)r are basically equal, indicating 
that the shorter length precludes reconnection. These results are consistent with the results of (HO, where Wg and Wt 
have the same respective meanings as in this paper. In that earlier work, it was found that the growth rates 7 scale as the 
appropriate fractional power of T] for Case I and scale linearly with r] for Case II. Both Q and A(^/ showed localized but 
nonsingular structures in the areas in which the corresponding periodic system has hyperbolic closed field lines or (for 
larger amplitude modes) chaotic behavior. On the other hand, the squashing factor shows similar behavior for Cases 
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I and 11. In particular, it shows almost no variation near the mode islands for small amplitude modes and localized 
behavior for larger mode amphtudes, with no clear signature of the different mode structure in Case I and Case 11. 



Acknowledgments 

We wish to thank V. Titov, E. Zweibel, and V. Mirnov for valuable suggestions. This research was supported by the 
DOE Office of Science, Fusion Energy Sciences and performed under the auspices of the NNSA of the U.S. DOE by 
LANL, operated by LANS LLC under Contract No. DE-AC52-06NA25396. 



Appendix: Issues related to the squashing factor Q 

As discussed in Sec.|2] the squashing factor is found in the following manner. First, let / be the Jacobian matrix 
for the map x ^ X which takes field lines from one surface, where Bn < 0, to the other surface, where Bn > 0. These 
surfaces are z = —L and z = L, respectively, in our example. For the special case B^q = const., we have D = det(/) = 1. 
The squashing factor is then simply the trace of J^J, as given by Eq. ^ with D = I. For the more general flux 
preserving case with B^ not constant, the determinant is D = ,X2, —L) ,X2,L) 9^ 1. It is tempting to change 
to canonical variables, i.e. variables u and U in which the equations of motion are canonical and the map is area 
preserving. This is always possible by the Darboux theorem, and it may appear at first blush that the situation is 
simplified because the determinant in these variables is unity. 

However, for the physical problems we consider, there is a natural metric on the two surfaces in the original 
variables, obtained from the Euclidean metric in 3D, so lengths are written as (5x, G5x). (For our example, the surfaces 
z = —L and z = LsLve flat, so this metric is the identity, G = I.) For general G, the relevant Rayleigh quotient is 

(5x,G5x) ' ^^^^ 

where Jij = dXi/dxj. The relevant eigenvalues are the eigenvalues of G~^J^GJ. 

The metric in the new canonical variables inherits the metric (5u, G5u) = (5x, G5x). This implies that G satisfies 
M^GM = G, where M is the tangent map from 5x to 5u (and from 5X to 5U). The new Rayleigh quotient is 

_ {5u,FGJ5n) 

(5«,G5«) ' ^^^^ 

which equals Qr. Therefore, if we transform to canonical variables and inherit the metric from the (physically 
motivated) metric on the original surface, we do not gain any simplicity from the relation D=\. 

This issue is related to the issue that arises for finite time Lyapunov exponents (FTLE), where the interval — L <z<L 
is replaced by the time interval < ^ < T, and x^X{t = T). If we change variables but continue to use the Euclidean 
metric, the FTLE 



h{T) = iln 



(5X(T),5X(T)) 



(26) 



(5x, 5x) 

are not coordinate invariant. The infinite time (largest) Lyapunov exponent, defined as /z = limr ^00 h{T), is coordinate 
invariant. The more natural question to ask is whether there is a natural metric. If so, then the metric changes upon a 
change of variables x ^ u and X ^ U so the FTLE are invariant as above. If there is no natural metric, then the FTLE 
depend on the metric, i.e. the Rayleigh quotient (5x,/^G/5x) / (5x, G5x) depends on the metric G. The infinite time 
Lyapunov exponents, however, are independent of G. 

A final issue relates to the angle between the magnetic field and the two surfaces. If B is far from normal to the 
surface, a circular flux tube projects to an elongated ellipse on the surface, and this is not related to reconnection. This 
issue was dealt with in Ref. 1 18|. For the model used in this paper, B^o ^ ^00 and the magnetic field is almost normal 
to the surfaces z = —L and z = L, so that this issue is not important. 
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